Mathematical properties of a semi-classical signal 
analysis method: noisy signal case 
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Abstract — Recently, a new signal analysis method based on a 
semi-classical approach has been proposed 1 1 1. The main idea in 
this method is to interpret a signal as a potential of a Schrodinger 
operator and then to use the discrete spectrum of this operator 
to analyze the signal. In this paper, we are interested in a 
mathematical analysis of this method in discrete case considering 
noisy signals. 

I. INTRODUCTION 

Recently, a new signal analysis method based on a semi- 
classical approach has been proposed [1|. We refer to this 
method SCSA for Semi-Classical Signal Analysis. The main 
idea in the SCSA is to interpret a signal as a potential of a 
Schrodinger operator depending on a semi-classical parameter 
0. It is well-known that if the potential is in the Faddeev 
class [3 1, then it can be expressed using a sum of the squared 
eigenfunctions associated to the negative eigenvalues charac- 
terizing the discrete spectrum of the Schrodinger operator and 
an integral involving the continuous spectrum. Similarly to 
the other standard approximation methods, by truncating the 
expression of the potential, the sum part is taken as an estimate 
of the potential. The proposed estimate depends on the semi- 
classical parameter. It has been shown that by reducing this 
parameter the estimation of the signal by the SCSA can be 
improved (T). 

Promising results have been obtained when applying the 
SCSA to arterial blood pressure. More than a satisfactory 
estimation of the pressure signals, this method introduced new 
spectral parameters that seem to contain important physiolog- 
ical information |4|, |5|. Moreover, a recent study has shown 
that the SCSA parameters could be useful in the estimation of 
some physical parameters related to turbomachinery features 
0. It has been also confirmed by some tests that the SCSA 
method is robust with respect to corrupting noises. Hence, 
it can be considered as a filter. However, the mathematical 
analysis of the SCSA in discrete noisy case has not been 
considered yet, comparing to other signal analysis methods 
(see, e.g., Q, 0, ®, (H, 031). 

We propose in this paper to study the mathematical prop- 
erties of the SCSA in discrete case considering noisy signals. 
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We propose also to study how to choose an appropriate semi- 
classical parameter to analyze a signal. In Section|IIl we recall 
the methodology and some mathematical properties of the 
SCSA method in the continuous case. Then, we study this 
method in discrete case in Section Hill In Section IIVI the 
SCSA method is studied in discrete noisy case. An a-posteriori 
error bound of the noise error contribution is given. In Section 
rvl we show how to choose an appropriate parameter for the 
SCSA method in discrete noisy case. Moreover, we can see 
its efficiency and stability with corrupting noises. Finally, we 
give some conclusions and perspectives for our future work in 
Section [Vj 

II. PRELIMINARY 

In this section, we are going to recall the methodology and 
some mathematical properties of the SCSA Q). 

A. Methodology 

Let us consider the following Schrodinger operator 

Hh{y):= - h2 £t- y > (i) 

where h £ is a semi-classical parameter [2], y £ C\(M.) := 
|V| \V(x)\ (1 + |x|) dx < ooj, and £{{R) is called the 



Faddeev class yj. Then, it is well known that y can be 
reconstructed as follows IPT21 



y(x) 



n=l 

lim \ f (*k 

b^+oo b J Q \tt 



4/i K nh ipl h (x)- 



kR r (i)h(k) f ±h {k,x)dk\ da, a. 



(2) 



where —K„. are the negative eigenvalues of Hh(y) with 
Ki h > K2 h > ■■ ■ > kn h > 0, Nh denotes the number of 
the negative eigenvalues, and %jj nh £ Ti. 2 {R) (H 2 (R) being the 
Sobolev space of order 2) are the associated C 2 -normalized 
eigenfunctions such that 



H h (y)i/j„ 



(3) 



Moreover, R r n)h(') is the reflection coefficient and fh±(-, ■) 
denote the Jost solutions defined as the unique solutions of 
the Schrodinger integral equation at ±00 respectively. 

Let us mention that the eigenpairs {n nh ,ip nh } for rih — 
1, • ■ • , Nh can be numerically calculated, but R r (i)h and fh± 
usually can not be calculated. Hence, similarly to the other 
standard approximation methods, by truncating (f2]) we propose 
to take the sum part as an estimate of y. Moreover, since this 

sum part Ah n nh i/% h is a positive definite function on R, 

n=l 

it is necessary to assume y to be positive definite too. Then, 
we give the following proposition. 

Proposition 1: JT] Let y be a positive definite function 
belonging to £}(R), then it can be estimated by 



Vh{x) 



N h 

:=4h^2K nh tpl h (x), 



(4) 



where h > 0, — are the negative eigenvalues of the 
Schrodinger operator Hh(y) defined by (Q~|i with K\ h > H2 h > 
■ ■ ■ > kn h > 0, Nh denotes the number of the negative 
eigenvalues, and ip nh are the associated C 2 -normalized eigen- 
functions. □ 
Consequently, the proposed estimation only depends on the 
parameter h JT]. 

B. Some properties of the SCSA method 

Now, it is natural to consider the convergence of this 
proposed method which is shown in the following proposition. 

Proposition 2: JH Let y be a real valued function satisfying 
the following condition 



yeB:= 
V eC\(M)\ Vi€l, y{x) >0, 



d m y 
g x m 



E C\R), to = 1,2 
(5) 



Then, we have lim^o \\yh - vWc^r) = °- 

Moreover, the number Nh of the negative eigenvalues of 
Hh(y) is a decreasing function of h. □ 

According to the previous proposition, we can improve our 
estimation by reducing the value of h. Moreover, it was shown 
in HI that if there exists an h such that is a reflectionless 
potential of the Schrodinger operator Hh(y), then yh is an 
exact representation of y. 

Let us recall that the study of the Schrodinger operator in the 
case where h — > is referred to the semi-classical analysis 0. 
Consequently, we call the proposed signal estimation method 
Semi-Classical Signal Analysis ( SCSA ) method. 

III. The SCSA method in discrete case 

In order to apply the SCSA method, we need to solve 
numerically the following Schrodinger eigenvalue problem 



H h (y)ip h = \ h iph, 4>h € H 2 (M.). 



Hence, we will study the SCSA method in discrete case in 
this section. 

We assume that y d is a discrete signal defined on an interval 
I = [a, b] C K, which can be considered as the restriction of 
a function y satisfying the condition given in ([5), such that 



(7) 



where Xj = a + (j — 1) Air for j = 1, • • • , M, are equidistant 
points with the distance between two consecutive points Ax = 
We denote y d } as the value y d {xj) for j = 1, • • • , M. 
The discretization of the Schrodinger eigenvalue problem 
given in © leads to the following eigenvalue matrix problem 



Xh^h 



where tyy 



*l>h,l,lph; 



Ad 



■ ■ ■ ,1ph,M-l^h,M 

-h 2 D 2 -dmg(Y d ), 



(8) 



(9) 



diag(F d ) is a diagonal matrix whose elements are y d for j = 
1 , ■ ■ • , M, and D2 is a second order differentiation matrix 
given by a discretization method for differential equations [ 13 1, 
lfT4l . which is independent of h. 

Let us consider the negative eigenvalues — ZsjL of Af with 
Ki h > ■■■ > Kjv h > 0' where Nh is the number of the 
negative eigenvalues of A4 with < Nh < M. We denote 



their associated eigenvectors by $ nh 
for nh = 1, • • • , Nh- Moreover, we assume that 

M 



i.e. 



1 



(10) 



Thus, according to Proposition Q] we can construct an estima- 
tion of y by the SCSA method in this discrete case as follows 



n=l 



for j = 1, ■ ■ ■ , M. By writing the following equality 



y 



t = {y d - vh) + (yn - y d h ) , 



(ii) 



(12) 



(6) 



we can see that the estimation error in this discrete case for 
yd can be divided into two sources: 

1) the estimation error for y^. y d — yh, which corresponds 
to the truncated integral part in (0, 

2) the discrete numerical error for y d : yh — y d , which 
is produced by the discretization of the Schrodinger 
eigenvalue problem. 

It is shown in Proposition [2] that the number of the negative 
eigenvalues of the Schrodinger operator is decreasing with 
respect to h. This property is generalized to the discrete case 
in the following proposition. 

Proposition 3: We assume that the matrix A d — —h 2 D2 — 
diag(y d ) defined in (O satisfies the following conditions: 

(CI): Z?2 is symmetric, 

(C2): Z?2 is negative definite, 



(C3): the number of zeros in the diagonal of diag(F d ) is 
equal to h with < h < M — 1. 
Let us denote the number of the negative eigenvalues of 
by Nh, then we have 



V/i > 0, 0<N h <M-n, 



(13) 



where A/ is the size of the matrix At. Moreover, we have 




N h = M 



for < h < 



di ' 



(14) 



(15) 



where y d +1 (resp. yfj) is the smallest (resp. largest) strictly 
positive element in the diagonal of diag(F d ), and d\ (resp. 
g?m) is the largest (resp. smallest) eigenvalue of □ 
Proof. Since the matrices D 2 and diag(y d ) defined in (0 both 
are symmetric, A^ is also symmetric. The eigenvalues of At, 
—D2 and — diag(F d ) can be denoted and ordered respectively 
as follows 



-vt 



Aft,M < • • ■ < Xh,u d M <■■■ < di, -y M < ■■■ < 
Then, according to Weyl's theorem (see |fT31 ) we get 

h 2 d M - yf < Xhj < h 2 d 1 - yf, (16) 

for j = 1, • • ■ , M, If n = 0, then ( fT3l can be directly obtained. 
By using (C2) and (C3) we obtain 

V/i > 0, < h 2 d M < A/ij, (17) 

for j = 1, • • ■ , n, and 



V^>0, Xhj <h 2 d 1 -y a J 



(18) 



for j = n+ 1, • • • , M. Hence, ( fT3l can be deduced from ( fTTI i. 
Moreover, if < h < M — 1, then by using ( TT8l we obtain 



that for any h < y ^ti, A/jj < 0, for j = n + 1, • • ■ , M. It 
yields Nh > M — h. Then, (fl4l can be obtained by using ( fT3] >. 
Finally, this proof can be completed by solving the following 
inequality 

0<h 2 d M -y d M <X h ,M- (19) 

□ 

Corollary 1: We assume that the matrix A^ defined in (0 
satisfies the conditions (C1)-(C3) given in Proposition|3] Then, 
V h > 0, 3 h! with < h! < h, such that 



N h , > Nh, 



(20) 



where N^ and Nh denote the number of the negative eigen- 
values of Af, and Af respectively. □ 
Proof. Let h > 0, by using (fT3T > we have < N h < M - h. 
Then, by using ( Tl3T > we get 



V < h! < min 



Thus, this proof is completed. 




, N h - > N h . 



(21) 



□ 



IV. ERROR ANALYSIS IN NOISY CASE 

In this section, we are going to consider the SCSA method 
in discrete noisy case. Moreover, an a-posteriori error bound 
of the noise error contribution for the SCSA method will be 
given. 

A. The SCSA method in discrete noisy case 
From now on, we assume that 

4 1 „ 



V 



y 



(22) 



is a noisy observation of the discrete signal y d defined in 
(Q, where the corrupting noise w is an identically distributed 
sequence of random variables with an expected value [i and 
a variance a 1 (a g R+). We denote yj 7 and Wj as the values 
y ru (xj) and znixj) respectively, for j = 1, - ■ ■ , M. 

In order to apply the SCSA method in this discrete noisy 
case, by substituting A^ in I© by A™ we need to solve the 
following eigenvalue matrix problem 



where 
diag^ 

i = V 



Al = -h 2 D 2 - diag(F ro ), 

is a diagonal matrix whose elements are yj 
, M. Hence, according to © we obtain 



At. 



Ad 



diag(W), 



(23) 

(24) 
7 for 

(25) 



diag(VF) is a diagonal matrix whose elements are zuj for j = 
1,-.. ,M. 

Let us denote the negative eigenvalues of A^ by — k 2 nh 



with ki 



> 



> h iw > for n h = 1, 



, Nh, where 



Nh is the number of the negative eigenvalues of A™. Similar 
to ( [Tol l, their associated eigenvectors are denoted by ^ nh = 

1 



A I 



with 



Thus, according 



to Proposition [T] an estimation of y can be given by the SCSA 
method in this discrete noisy case as follows 



Ah^k rih ^ 2 nh 

n=l 



for j = 1, 



Vh 



, M. By writing the following equality 

= (y d Vh) + (vh - y d h ) + (y d h - yf) , 



(26) 



(27) 



we can see that the total estimation error in this discrete noisy 
case for y^ can be divided into three parts: 

1) the truncated error for y^: y d — yh, 

2) the discrete numerical error for yt: yh — y d , 

3) the noise error contribution for y^\ y d — y™. 
Proposition [4] illustrates a property on the number of the 

negative eigenvalues in the discrete noisy case. 

Proposition 4: We assume that the matrix A\° = —h 2 D2 — 
diag(Y ra ) defined in (GUi satisfies the conditions (C1)-(C2) 
and the following condition 

(C4): the number of the positive elements in the diagonal 
of — diag(F ro ) is equal to h with < n < M — 1. 



Let us denote the number of the negative eigenvalues of 
by Nfi, then we have 



Vh>0, < N h < M-n, 



(28) 



for nh = 1, • • ■ ,Nh- By using the Bienayme-Chebyshev 
inequality we get that for any real number 7 > 0, 



Pr(H-A*| <7<0 > 1- -5, 



(34) 



A^ = 0, for h> J^, 

V «M 



N h =M-h, for < /i< 



di ! 



(29) 
(30) 



where (resp. y^) is the smallest (resp. largest) strictly 

positive element in the diagonal of diag(y ro ), and d\ (resp. 
<1m) is the largest (resp. smallest) eigenvalue of D^. 

Moreover, we have V h > 0, 3 h! with < hi < h, such 
that N h i > Nh, where Nh' and Nh denote the number of the 
negative eigenvalues of AJ^, and respectively. □ 
Proof. This proof can be completed in a similar way to the 
one of Proposition [3] □ 

B. Analysis of the noise error contribution 

In this subsection, we are going to study the noise error 
contribution in the SCSA method by providing an a-posteriori 
error bound in the following proposition. 

Proposition 5: Let y ro be a discrete noisy signal defined as 
in d22i i. y^ and yi be the estimations of y given by (fTTT i and 
(|2oT l respectively. Moreover, we assume that 

(CI): the matrix D2 given in (O and ( f24T > is symmetric, 

(C5): the numbers of the negative eigenvalues of and 
A™ are equal, i.e. N h = N h , 

(C6): kl h < 2 kl h , for n h = 1, • • • , N h . 
Then, an a-posteriori error bound for the noise error contribu- 
tion in y^ can be given as follows 



\Vh 



p 
< 



4h 

7Kc 



N h 

E 

n=l 



1R r 



B 1 



(31) 



where a = max ( 



|// — y<t|, I// + 70- 1) with /i and c 2 being 

the expected value and the variance of zuj respectively for 

p 

j = 1, • • • , M, and c < d means that the probability for c to 
be smaller than d is p with p = 1 — -4j an d 7 € R+. □ 
In order to prove this proposition, we need the following 
lemma. 

Lemma 1: By giving the conditions (CI), (C5), and (C6), 
we have 



v 
< 



B 1 
y2 K n , 



for nh = 1, • • • , iVft, where S^ CT = max (|/x - 7*7 1, |/i 



(32) 

7cr|). 
□ 



Proof. According to the condition (CI), the matrixes and 
A™ are both symmetric. Moreover, since the matrix diag(M / ) 
is diagonal, its eigenvalues are Wj for j = 1, • • • , M. Hence, 
by using Weyl's theorem (see [ 15 1) and (l25T l with the condition 
(C5), we obtain 



\ K n h 



< max \xuj 

l<j<M J 



(33) 



i.e. the probability for tUj to be within the interval 
]/x — 7a , /x + 7<r[ is higher than 1 — A. Consequently, by 
denoting CT = max — 70"], |/z + ~fcr\) and using ( f33b we 
obtain 

(35) 



for = ,Nh- According to the condition (C6) we 

obtain 



0<^ <min(^ )/i 2j 



(36) 



Then, the utilization of the mean value theorem gives us 



1 "/1 



" 1 1 " 2 



1 \ K n h K n h I 



(37) 



Then, this proof can be completed by using ( 1351 ). □ 
Proof of Proposition |5J By using the condition (C5), ( fTTT i and 
( l26l >. we get 

n=l 

for j = 1, • • • , M. Let us denote E nh = [e nh ,i, ■ • • , e„ ft , A f] T 
with 

(39) 



Then, we obtain 



(40) 



By calculating the norms 
be bounded as follows 



\E nh \\ 2 < 



and 



Hence, by using 



and (f4TT >. we get 



l-Sn h || a can 



(41) 



C-2/ft|| 2 <4ftEl|S„J 



4/? 



(42) 



v n— 1 



Finally, the proof can be completed by using Lemma Q] □ 
The convergence of the SCSA method in the continuous 
noise-free case is shown in Proposition^ where the estimation 
proposed by the SCSA method can be improved by reducing 
the value of h. Now, it is interesting to study the efficiency 
of this method in discrete noisy case. Especially, we need to 
know the influence of h on the noise error contribution in the 
SCSA method. A natural idea is to study the influence of the 



parameter h on the noise error bound given in Proposition [5] 
so as to deduce the one on the noise error contribution. On 
one hand, because of the term h this error bound seems to be 
increasing with respect to h. On the other hand, it is shown 
in Proposition @] that the number of the negative eigenvalues 
Nh can be decreasing with respect to h. Consequently, it is 
impossible to intuitively know the influence of h. However, it 
can be studied numerically as shown in the next section. 

V. NUMERICAL RESULTS 

In this section, by taking a numerical example we are going 
to show how to choose an appropriate value of h for the 
SCSA method in discrete noisy case. Moreover, we can see 
its efficiency and stability with corrupting noises. 

We assume that y ro is the discrete noisy observation of y 
defined in (f22t . where 



V d (xj) = y(xj) = sech 2 (xj 



(43) 



with xj e / = [0, 12] and Air = 1CT 2 . Hence, M is equal 
to 1201. Moreover, we assume that the noise za is simulated 
from a zero-mean white Gaussian iid sequence. The variance 
a 2 of w is adjusted in such a way that the signal-to-noise ratio 
is equal to lldB. We can see the original signal y and its noisy 
observation y™ in Figure Q] Let us recall that this sech-squared 
function is well known in the quantum physics theory as the 
Poschl-Teller potential of the Schrodinger operator [16|. 

In order to estimate y by using (126) . we propose to use 
a Fourier pseudo-spectral method |[T3l to solve numerically 
the Schrodinger eigenvalue problem defined in ©. Thus, the 
second order differentiation matrix D 2 is given as follows [ 14 1: 
If M is even, then 



D 2 (k,j) 



If M is odd, then 

D 2 (k,j) 




3A 2 



J_ 

12 ' 



(Ax) 



■ e _j 1 cot(i 



for k = j, 
for k j. 



for k = j, 
for k 7^ j. 



with A = Let us mention that D 2 is symmetric and 
negative definite. Then, we use the Matlab routine eig to solve 
the eigenvalue matrix problems defined in ([8]) and < f23b - 

It is shown in Subsection |IV-B| that the total estimation error 
for y^ comes from three parts. However, since the estimation 
yh can not be calculated in the discrete case, we only consider 
the estimation error in the discrete noise-free case and the 
noise error contribution 



y d - C = (v d - vt) + (y d h - vu) 



(44) 



In order to see the influence of h on the total estimation 
error, we show the variations of \\y d — y^\\ 2 in Figure [3(a)] 
which is represented by the black solid line. We can see that 
\\y d ~ y^\\ 2 nas a mm i mum at h = 0.4 and a local minimum 
at h = 0.7. Thus, we can take the optimal value h — 0.4 
for y^ so as to produce a minimal total estimation error. The 



estimation obtained by using with h = 0.4 is given in 
Figure [2] Consequently, we can see that the SCSA method is 
accurate and robust with a corrupting noise. Hence, it can be 
considered as a filter for noisy signals without delays. 

The previous analysis is based on the knowledge of y d 
which is usually unknown in the practice work. If we use 
— y™\\ 2 to study the influence of h on the total estimation 
error, then we generally can not find an optimal value of h. 
This can be explained by the green solid line in Figure |3(a)| 
which corresponds to the different values of ||j/ ro — y^|| 2 - I n 
order to solve this problem, we propose to use a second-order 
Butterworth filter ifTTIl which is given as follows 



H(s) 



s 2 + 2w c s + w 2 . 



(45) 



where the cutoff frequency w c is set to w c = 0.01. 

H(-) is a classical low -pass filter which can be used to 
attenuate the corrupting Gaussian noise in The filtered 
signal is represented by the black dotted line in Figure [2] Since 
this filter produces a delay to y™ , we also apply it to y™ . 
The filtered signals are denoted by and fy^ respectively. 
Then, we use the influence of h on — fVh\\ 2 to deduce 



the one on 



1 2/" — 2//? 1 1 2 wn i cn corresponds to the estimation 
error in the discrete noise-free case. We can see in Figure 



fVh 



and 



|3(a)| the relation between the variations of ||/y 
\\y d — y d \\ 2 which are represented by the red dotted line and 
the blue dash-dotted line respectively. Hence, we can observe 
that they have the same variation with the same local minimum 
and the same local maximum. 

Now, we study the influence of h on the noise error 
contribution. The variations of ||y^ — y^ || 2 is shown in Figure 
|3(b)| Moreover, we can verify that the number of the negative 
eigenvalues Nh and Nh are equal. We can see their variation 
with respect to h in Figure 3(d)| Consequently, according 
to Figure |3(b)| and Figure |3(d) we can observe that when 
Nh is equal to 1 the noise error contribution is increasing 
with respect to h. However, when the value of Nh increases, 
the noise error contribution is decreasing with respect to h. 
Consequently, a small value of h can produce a large noise 
error contribution. 

We are going to use Proposition [5] to deduce the variation of 
the noise error contribution with respect to h. Since the noise 
is assumed to be a zero-mean white Gaussian iid sequence, 
by using the well known three-sigma rule we obtain that 



99.7% 

max kx7, < 3<7. 

l<j<M J 



(46) 



The variation of the noise error bounds given in Proposition 
is shown in Figure [3(c)| Although this error bound is not sharp, 
its variation is similar to the one of the noise error contribution 
shown in Figure |3(b)| where there are local minimums at h = 
0.6 and h — 0.4, and a local maximum at h = 0.5. 

Finally, by combined with the variation of — fVh\\ 2 

shown in Figure [3(a)] we can choose ft, = 0.4or/i = 0.7in 
our estimation so as to minimize the total estimation error. 



Original signal and its noisy observation 

1.4r " , 




Time 

Fig. 1. The discrete signal y d and the discrete noisy signal with SNR = 
lldB. 



Original signal and estimated signals 

1.2 r " i 




2 4 6 8 10 12 

Time 

Fig. 2. The discrete signal y d , the signal jy^ filtered by Butterworth filter, 
and the estimate y^ obtained by the SCSA method with h = 0.4. 




(a) The I 2 norms of different errors. (b) The noise error contribution. 
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(c) The noise error bound (d) The number of negatives eigenval- 

ues. 

Fig. 3. Results of different variations obtained for h = 0.2, 0.3, ■ • * , 1.9, 2. 



VI. CONCLUSION 

In this paper, the SCSA method recently introduced for 
signal analysis is studied in discrete noisy case. Some mathe- 
matical properties of the negative eigenvalues of a Schrodinger 
operator are given in discrete noise-free case and discrete noisy 
case respectively. An a-posteriori error bound of the noise 
error contribution is given which is based on the expected 
value and the variance of a corrupting noise. By taking 
a numerical example, we show the influence of the semi- 
classical parameter on different sources of errors in the SCSA 
method. Moreover, it is shown that the SCSA method is 
accurate and robust against corrupting noises. Hence, it can 
be considered as a filter without involved delays. Finally, we 
study how to choose an appropriate semi-classical parameter 
without knowing the original signal. The comparison to other 
signal analysis methods like Fourier transform or the wavelets 
will be done in a future work. Moreover, the SCSA method 
will be extended for time derivatives filtering, which is still 
an open problem, such that this method can be useful in more 
applications in signal processing and automatic control. 

References 

[1] T. M. Laleg-Kirati, E. Crepeau and M. Sorine, Semi-classical signal 

analysis, to appear in Mathematics of control signals and systems 

(MCSS) journal, 2012 
[2] M. Dimassi and J. Sjostrand, Spectral asymptotics in the semi-classical 

limit, Cambridge U Press, 1999. 
[3] L. D. Faddeev, Properties of the S-matrix of the one-dimensional 

Schrodinger equation, Trudy Mat. Inst. Steklov, 73, pp. 314-336, 1964. 
[4] T. M. Laleg-Kirati, C. Medigue, F. Cottin and M. Sorine, Arterial blood 

pressure analysis based on scattering transform II, in Proc. EMBC, Lyon, 

France, August 2007. 
[5] T. M. Laleg-Kirati, C. Medigue, Y. Papelier, F. Cottin and A. Van de 

Louw, Validation of a semi-classical signal analysis method for stroke 

volume variation assessment: A comparison with the PiCCO technique, 

Annals of biomedical engineering, vol. 38, no. 12, pp. 3618-3629, 2010. 
[6] F. Eleiwi, T. M. Laleg-Kirati, S. Khelladi and F. Baki, A Semi-Classical 

Signal Analysis Method for the Analysis of Turbomachinery Flow 

Unsteadiness, World Academy of Science, Engineering and Technology, 

vol. 59, pp. 215-218, 2011. 
[7] M. Fliess, Analyse non standard du bruit, C.R. Acad. Sci. Paris Ser. I, 

342, pp. 797-802, 2006. 
[8] D. Y. Liu, O. Gibaru, W. Perruquetti, M. Fliess and M. Mboup, An 

error analysis in the algebraic estimation of a noisy sinusoidal signal. In: 

16th Mediterranean conference on Control and automation (MED'08), 

Ajaccio, France, 2008. 
[9] D. Y. Liu, O. Gibaru and W. Perruquetti, Error analysis for a class 

of numerical differentiator: application to state observation, 48th IEEE 

Conference on Decision and Control, Shanghai, China, 2009. 
[10] D. Y. Liu, O. Gibaru and W. Perruquetti, Error analysis of Jacobi 

derivative estimators for noisy signals. Numerical Algorithms, 58, 1, 

pp. 53-83, 2011. 

[11] D. Y. Liu, O. Gibaru and W. Perruquetti, Convergence Rate of the Causal 

Jacobi Derivative Estimator. Curves and Surfaces 2011, LNCS 6920 

proceedings, pp. 45-55, 2011. 
[12] D. J. Kaup and A. C. Newell, Solitons as particles, oscillators, and in 

slowly changing media: a singular perturbation theory, Proc. R. Soc. 

Lond. A. 361, pp. 413-446, 1978. 
[13] J. P. Boyd, Chebyshev and Fourier Spectral Methods, Dover Publication, 

Inc., 2nd edition, 2000. 
[14] L. N. Trefethen, Spectral Methods in Matlab, SIAM, 2000. 
[15] G. Stewart and J. G.Sun, Matrix Perturbation Theory, Academic Press, 

New York, 1990. 

[16] G. Olafsson and S. Zheng. Function spaces associated with Schrodinger 
operators: the Poschl-Teller potential. Journal of Fourier Analysis and 
Applications, 12(6), pp. 653-674, 2006. 



G. Bianchi and R. Sorrentino, Electronic Filter Simulation & Design. 
New York: Mc Graw-Hill, 2007. 



